Meinardus' theorem on weighted partitions: 
extensions and a probabilistic proof 



Boris L. Granovsky 

Department of Mathematics, Technion-Israel Institute of Technology, Haifa, 

32000, Israel 

Dudley Stark 

School of Mathematical Sciences, Queen Mary, University of London, London El 

4NS, United Kingdom 

Michael Erlihson 

Department of Mathematics, Technion-Israel Institute of Technology, Haifa, 

32000, Israel 



Abstract 

The number c n of weighted partitions of an integer n, with parameters (weights) 
bk, k > 1, is given by the generating function relationship Y^=Q c n zU = IlfcLiU ~~ 
z k )~ bk . Meinardus(1954) established his famous asymptotic formula for c n , as n — > 
oo, under three conditions on power and Dirichlet generating functions for the se- 
quence bk- We give a probabilistic proof of Meinardus' theorem with weakened third 
condition and extend the resulting version of the theorem from weighted partitions 
to other two classic types of decomposable combinatorial structures, which are called 
assemblies and selections. 
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1 Summary 



In this paper, we combine Meinardus' approach for deriving the asymptotic 
formula for the number of weighted partitions with the probabilistic method 
of Khintchine to develop a unified method of asymptotic enumeration of three 
basic types of decomposable combinatorial structures: multisets, selections 
and assemblies. As a byproduct of our approach we weaken one of the three 
Meinardus conditions. In accordance with these two objectives, the structure of 
the paper is as follows. Section 2 presents Meinardus' asymptotic formula, the 
presentation being accompanied by remarks clarifying the context of the three 
conditions of Meinardus' theorem. In section 3 we state our main result which 
are asymptotic formulae for numbers of multisets, selections and assemblies. 
Sections 4, 5 and 6 are devoted to the proof of the main theorem, including 
the unified representation of basic decomposable random structures, which is 
the core of the probabilistic method considered. In Section 7 we discuss the 
striking similarity between the derived asymptotic formulae. 



2 Meinardus' theorem 



The Euler type generating function /W for the numbers c£\ n > 1 of weighted 
partitions of an integer n, with parameters bk > 0, k > 1 is 

oo oo 

/ (1) W:=E t 4 1)2B = II( 1 -«*)A W < L C 1 ) 

n=0 k=l 

In this setting, bk is interpreted as a number of types of summands of size 
k. (For example, one can imagine that coins of a value k are distinguished 
by bk years of their production). It is also assumed that in a partition, each 
summand of size k belongs to one of the bk types. In the case bk = 1 for all 
k > 1, is the number of standard (no n- weighted) partitions of n (with 
Co = 1), while the case bk = k, k > 1 conforms to planar partitions, studied 
by Wright, see pQ and the recent paper [T7| by Mutafchiev. The study of the 
asymptotics of the general generating function ([T]) was apparently initiated 
by Brigham who obtained in [6] the asymptotic formula, as n — > oo for the 
logarithm of the function, using the Hardy- Ramanuj an asymptotic technique. 
Meinardus' approach ([16]) to the asymptotics of is based on considering 
two generating series for the sequence bk > 0, k > 1: the Dirichlet series D(s) 
and the power series G(z), defined by 

oo 

£(s) = E & ^~ S > s = a + it, (2) 

k=l 
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G(z) = J2hz k , \z\<l. (3) 

k=l 

We note that the function f^'(z) converges at the point \z\ < 1 if and only if 
the same is true for the function G(z) (see e.g. Lemma 1.15 in [7]). 

Meinardus ([IS]) established the following seminal asymptotic formula for cffl, 
which is presented in pQ. We denote by and the real and imaginary 
parts of a number. 

Theorem 1 (Meinardus) Suppose that the parameters bk > 0, k > 1 of 
weighted partitions meet the following three conditions: 

(i) The Dirichlet series |J|) converges in the half-plane a > r > and there 
is a constant < Co < 1, such that the function D(s), s = a + it, has an 
analytic continuation to the half-plane 

H = {s:a> -C } (4) 

on which it is analytic except for a simple pole at s = r with residue A. 

(ii) There is a constant C\ > such that 

D(s) = O (\t\ Cl ) , t-^oo (5) 



uniformly in o > —Co- 
(Hi) There are constants C<i > 0, e > such that the function 

g{r) := G(exp(—r)), r = 5 + 2iria, 5 > 0, a e R, (6) 

satisfies 

®(9(t)) - g(5) < -C 2 5~% |argr|>^, ^ \a\ < 1/2, (7) 
for 5 > small enough. 
Then, as n — ► oo, 

c« ~ C^n^ exp Cn r /( r+1 )(l + -) (AT(r + l)((r + l)' 1 ' 



where 

and 
where 



2D(0) - 2-r 

K l = : ; 

2(1 + r) 

C (i) = e D'(o) ( 2vr (i + r ))-V2 (AT(r + l)C(r + 1)) 

1 - 2D(0) 
2(1 + r) 



K-2 
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Meinardus also gave a bound on the rate of convergence which we have omitted 
in the statement of Theorem [TJ 



At this point we wish to make a few clarifying comments on the three Meinar- 
dus conditions (i), (ii), (Hi). 

• The Ikehara- Wiener Tauberian theorem on Dirichlet series cited below tells 
us that condition (i) implies a bound on the rate of growth, as k — > oo of 
the coefficients bk of the Dirichlet series D(s) in fl2j). 
Theorem 2 (Wiener-Ikehara) (see Theorem 2.2, p. 122 in Tffl$ ) 

Suppose that the Dirichlet series D(s) = J2T=i a kk~ s is such that the 
function D(s) — -^y has an analytic continuation to the closed half-plane 
9?(s) > 1. Then, ' 



k=l 



au ~ An, n — ► oo. (9) 
We will use the fact that (jSJ) implies 

cifc = o{k), k — > oo. (10) 
To prove this, we rewrite as 



71-1 



1 n 1 

2J a k = —a n + - a k = A + e„, e n — > 0, n ^ oo, 



n ~i " " fc=i 



which gives 



-a n H (-4 + e n _i) = A + e n . 

n n 



Consequently, lim n ^ 00 (l/n)a n = 0. 

Now set ak = k~ r+1 b k , k > 1, where bk, k > 1 satisfy Meinardus' condi- 
tions (i) and (ii). Since C ,r > 0, the sequence a k obeys the conditions of 
the Wiener-Ikehara theorem, so that we get from (flO]) the bound: 

b k = o{k r ), k^oo. (11) 

Functions satisfying Meinardus' condition (ii) are called of finite order in 
the corresponding domain. It is known (see e.g. [20], p. 298) that the sum 
D of a Dirichlet series is a function of a finite order in the half-plane of the 
convergence of the series. Thus, condition (ii) requires that the same holds 
also for the analytic continuation of D in the domain Ti,. 
We show below that the condition (Hi) is associated with bounding the so- 
called zeta sum known from the theory of the Riemann Zeta function. In 
fact, 



$t(g(r)) -g(5) = -2 he~ kS sin 2 (nka), 5 > 0, aG 

k=l 
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which allows us to reformulate © as 



S 



2 V b k e~ kS sin 2 (7rfca) > C 2 5~ e , < — < \a\ < 1/2, (12) 

for S > small enough and some e > 0. 

The verification of condition (Hi) in the forthcoming Lemma [1] relies on 
the lower bound (jHj) below, for the sum I^ =1 sin 2 (7r/ca), aeR. This bound 
can be derived from the following bound on the zeta sum in the left hand 
side of (USD (see p2], p. 112, Lemma 1): 



p 

k=l 



< min < P. 



' 2 II a 



P > 1, a G 



(13) 



where || a || denotes the distance from a to the nearest integer. It follows 
from (pi]) that for all a E R 



2 shi 2 (irka) > P 
k=i 



p 

^ ' ^2nika 
k=l 



> P — min < P, 



2 a 



which is convenient to rewrite as 

2^sin 2 (7rAxt) > P |l -minjl, 

Under the assumptions in Meinardus' condition (Hi). 



a 



(14) 



/II a ||= \a\ < 1/2, |a|5 _1 > — . 

27T 



Setting in flU 



P = P(a,5) 



1 + 


air 1 


2 







> 1, 



(15) 



(16) 



where [x] denotes the integer part of x and 5 > is small enough, we get 
the desired bound, 



p g-i 
2]Tsin 2 (7rAxO > — , 

k=l 



(17) 



provided ( TTBl) holds. It follows from the above that for any fixed ko > 1, 
and any < e = e(8] k ) < 1/2, 



k=ko 



2 sm 2 (7rka) > f- - e) cT 1 := aT 1 , 



if 5 > is small enough and (TT5]) holds 
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In the proof of Lemma [T] below we will also use the fact that under the 
condition (fl5|) . the choice (fT6l) of P provides, 

It seems not to have been noticed that Meinardus' condition (Hi) is rather 
easily satisfied, as is shown in the following lemma. 

Lemma 1 Let the sequence {bk} be such that bk > pk r ~ l , k > k for some 
k > 1 and some constants p, r > 0. Then HUty is satisfied. 

Proof Because of f fl6|) . P > and therefore P > k for 5 > small 

enough. We have, 

oo p 

2 he~ k5 sin 2 (nka)> 2 ]T pk^e^ 5 sm 2 (nka) 

k=l k=ko 

P 

> 2pe- ps k r ~ l sin 2 (Tika) := Q. 

k=ko 

In order to get the needed lower bound on Q implied by (Tl2i) . we need to 
distinguish between the following two cases: (i) < r < 1 and (ii) r > 1. 
Applying ffl8l) and ffl9l) we have in case (z), 

Q > pe-^P^cS- 1 > pe~ d (P5) r - 1 c5- r > pe~ d d r ' 1 c5- r 

and in case (ii), 

Q > pe~ d k r ~ l cS-\ 

Therefore, ffl2l) is satisfied with e = r in case (i) and with e = 1 in case (ii). ■ 

We note that in [T7j the validity of condition (Hi) was verified in the particular 
case of planar partitions (bk — k, k > 1), via a complicated analysis of the 
power series expansion of the function 3t,(g(r)) — g(5). 

Example 1 Let bk = pk r ~ x , p, r > 0, k > 1. Such weighted partitions are 
associated with the generalized Bose-Einstein model of ideal gas (see [2TJ ) . In 
this case, D(s) = p((s — r + 1), where ( is the Riemann zeta function. Thus, 
D(s) has only one simple pole at s = r > with the residue A = p and it has 
a meromorphic analytic continuation to the whole complex plane C. These 
facts together with Lemma [1] show that all three of Meinardus' conditions 
(i), (ii), (ii) hold. In the case considered the values D(0) = p((l — r) and 
-D'(O) = pC'(l — r ) i n the asymptotic formula (IE]) can be found explicitly from 
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the functional relation for the function (, as is explained in [T7]. In particular, 
for standard partitions (p = r = 1), 

D(0) = C(0) = ~ D'(0) = C'(0) = -\ log2vr, 



while for planar partitions (p — l,r = 2), 



I>(0) = C(-l) = ~ , J D'(0)=C / (-l) = 2^ 



00 wlogw 

aw. 



g27TUI J 

For an arbitrary r > 0, the expressions for -D(O), -D'(O) include the integral 

u/'" 1 log w 



Example 2 The purpose of this example is to show that conditions (i) and 
(22) of Theorem [T] do not imply condition (Hi) in the same theorem. Let 



1, if4|Jfe 
0, if Aj(k. 



Let a = 1/4 in the sum X)&Li ^fc e 5fe sin 2 (7r/ca). Then, because for all either 
6 fc = or sin 2 (7rA;/4) = 0, 



^6 fc e-*sin 2 (7rA;/4) = 



k=l 



and therefore ffT2l is not satisfied. However, 



3=1 

which clearly satisfies the first two of Meinardus' conditions because the func- 
tion 4~ s is entire and |4 _s | = 4~ CT < 4 C ° for s G Tl, where 7i is given by 



3 Statement of the main result 



Our main result, Theorem [3] below, achieves two objectives: weakening the 
Meinardus condition (Hi) and extending the resulting version of the Meinardus 
theorem from weighted partitions to other two types of classic decomposable 
combinatorial structures. 
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We first recall that a decomposable structure is defined as a union of in- 
decomposable components of various sizes. It is known (see [2113] ) that the 
three types of decomposable combinatorial structures: multisets, which are 
also called weighted partitions, selections and assemblies, encompass the va- 
riety of classic combinatorial objects. Weighted partitions are defined as in 
the previous section, selections are defined as weighted partitions in which 
no component type appears more than once and assemblies are combinatorial 
objects composed of indecomposable components which are formed from la- 
belled elements. Each decomposable structure is essentially determined by the 
number of types of its indecomposable components having a given size k. We 
denote this number by bk for weighted partitions and selections and by m& for 
assemblies. In the case of assemblies we denote bk = rrik/k\, so that in all three 
cases h)~, k > 1 are parameters defining a structure. In what follows we will 
use the notation »W, i — 1,2,3 for quantities related to weighted partitions, 
selections and assemblies respectively. Given a sequence bk, k > 1, we define 
cW = for i — 1, 2 and define cffi = s^/n!, where s$ denotes in all three 
cases the number of combinatorial structures of type i having size n. 

Theorem 3 Suppose that the parameters bk, k > 1 meet Meinardus' condi- 
tions (i) and (ii) as well as the condition 

(Hi') For 5 > small enough and any e > 0, 

00 / r \ 

2 he- kS sin 2 (nka) > 1 + - + e | logo"), 
fc=i V 2 / 

V^<H<l/2, z = l,2,3, 
where the constants Af^ are defined by 

fkjs' */ i = 1 . 
M (i) = ^ 4, if i = 

U */ i = 3. 

T/ien t/ie asymptotics for c£\ i = 1,2,3, as n — > oo, are given respectively 
by Meinardus ' formula $8\), and by the formulae UZO) , WIS) below: 

c {2) ~ C {2) n~^ exp (n^Ti (l + I) (Ail - 2~ r )C(r + l)T(r + 1)) , (20) 
where 

C (2) = 2 D(p) ^ 27r (l + r )J ^r(r + l)(l -2- r )C(r + l)) 2 " +2 , 
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and 

cf ~ C (3) n"^ exp (n^rr (l + _) (AT(r + 1)) ^ , (21) 

/ \ _1 / 2 i 

C (3) = e D(0) 2tt(1 + r) (AT(r + l))w . 



Remark H-K Hwang ([H]) applied the approach of Meinardus to the study 
of the asymptotics of the number of summands, say u n , in weighted parti- 
tions and selections, which he called unrestricted and restricted partitions 
respectively. In the first local limit theorem for a properly scaled u n 

was obtained in [TT] under the three conditions of Meinardus. Regarding re- 
stricted partitions, the author claimed the same under Meinardus' conditions 
(i), (ii) and a condition similar to our (Hi'), but the proof contains an error 
in bounding the function Gg(r) on p. 109. 

Example 3 This example satisfies all three conditions of our Theorem [3X but 
does not satisfy condition (Hi) of Theorem CD Let bk, k > 1, be defined by 



12e 7 (^) if4/fc, 
12e 7 (50 + ^-2MM) if 4|fc and 16/A;, 

l 2e 7( 5 o + _ 2 iHS^! + if 16\k. 

Note that because < < e" 1 for x > 1, it follows that b k > 12e 7 (^) , for 
all k > 1. The Dirichlet series D(s), s = a + it, for this choice of converges 
absolutely for a > 1 and in this domain 

D(s) = 12e 7 (-(I- A~ s ) 2 C(s + 1) + 50 ■ ^((s)) , (22) 



where we have used the fact that ('(s + 1) = — J2kLi X ^^k~ s . It is well known 
that the function ((s) has a simple pole at s — 1 and that the Laurent expan- 
sion of ((s + 1) around s = is 

C(s + l) = -+7 + -"- , (23) 
s 

where 7 is Euler's constant. It follows from (1231) that the function £'(s + 1) 
has a unique pole at s = of order 2. As a result, we derive that in (1221 the 
first term in the parentheses is analytic in the whole complex plane C, while 
the function D in (12"2"1) is analytic in C except a simple pole at s = 1. It is also 
a known fact that the functions (, (' satisfy ([5]), from which we conclude that 
the same is true for the function D given by (|22|) . To show that condition (Hi 1 ) 
of Theorem [3] is satisfied, we note that, if 5 > and e > are small enough 
then 
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p 



J2 he~ kS sin 2 (nka) > J2 he~ kS sin 2 (rrka) 

k=l k=l 

>12eV d £^sin 2 (7rto) 
k=i k 



>12eV d ^£sin 2 (7rA;a) 

" k=3 

> 12eM !2i^) ( i_ er . (24) 

>61og(r 1 ), (2tt)- 1 (5 < |a| < 1/2, (25) 



where we have used ( [18]) and ([19]) at (1241) and the fact that 3.5 < d < 4 
in the last step. Since in the case considered r = 1, the condition (iW) is 
indeed satisfied for all three types of random structures. Finally, to show that 
condition (in) of Theorem [T] in the form f[T2"]) is not satisfied, we set a = 1/4 
in the left hand side of ffT2l) to obtain for 5 — > 0, 



f] 6 fc e" fc5 sin 2 (vr-) < 12e 7 £ ^fV^ ~ 12e 7 /°° ^e'^dx 

k=l ^ fc=l ^ Jl X 

i / — — -e x <ix + 12e^ / — — -e (ix = O log 2 (5 x ) 

h x Ji x v 



Example 4 Consider the assembly of forests, for which components consist of 
labelled linear trees. The number of such components on k vertices is rrik = k\ 
and so 6fc = 1, just as for integer partitions. The asymptotic number of labelled 
linear forests is thereby given by ( |2T[) in Theorem [3] with r — 1, A = 1. We 
note that the number of labelled linear forests on n vertices equals the number 
of path coverings of a complete graph on n vertices. 



4 A unified probabilistic representation for decomposable combi- 
natorial structures. 



It has been recently understood (see [T8I2T] ) that the three main types of 
decomposable random structures: assemblies, multisets and selections, are in- 
duced by a class of probability measures on the set of integer partitions, having 
a multiplicative form. Vershik ([21]) calls the measures multiplicative, while 
Pitman ([18]. [5]) refers to them as Gibbs partitions. Equivalently, in combi- 
natorics it is common to view the structures above as the ones generated by 
the conditioning relation (see [2]) or by the Kolchin generalized allocation 
scheme (|13j). Our asymptotic analysis is based on the unified Khintchine 
type probabilistic representation of the number of decomposable structures 
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of size n. Recall that we agree that the number of non labelled structures 
(weighted partitions and selections) is denoted by and respectively, 
and the number of labelled structures (assemblies) by n\c£> . In all three cases 
the probabilistic representation of c n is constructed as follows. Let / be a 
generating function of a sequence {c n } associated with some decomposable 
structure: 

f{z) = £ c n z n . 

n>l 

A specific feature of decomposable structures is that the generating function 
/ has the following multiplicative form: 

/ = n s k, 

k>l 

where Sj. is a generating function for some nonnegative sequence {dk(j), j > 
0, k > 1}, i.e. 

S k (z) = J2d k (j)z k \ k>l. (26) 

We now set z = e~ s+2nta , a G [0, 1] and use the orthogonality property of the 
functions e ~ 2man , n > 1 on the set a G [0, 1], to get 



Cn = e nS £ f (V 5+27riQ ) e- 2lTian da 

rl n / \ 
gn<5 / | | ( q ( a — S+2-Kia\ \ a —2nian 







II ( S k (e- 5+27Tia ) )e- 27Tian da, n > 1, 

k=l 



where 5 is a free parameter. We denote by 

n 

fn := II Sfc, n > 1, (27) 

k=l 

the truncated generating function. Next, we attribute a probabilistic meaning 
to the expression in the right hand side of (|27j) by defining the independent 
integer valued random variables Y k , k > 1 : 

d, (i)p~ Skj 

F(Y k =jk)= f , j>0, k>l. (28) 
and observing that 

is the characteristic function of the random variable 



Z n :=J2Yk- (30) 



k=i 
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We have arrived at the desired representation: 

c n = e n5 f n (e- 5 )F(Z n = n), n > 1. (31) 

In accordance with the principle of the probabilistic method considered, we 
will choose in fl3Tl) the free parameter 5 = 5 n to be the solution of the equation 

EZ n =n, n > 1, (32) 

after we show in the next section that for the three classic combinatorial 
structures the solution to (13"2l) exists and is unique. 



It can be easily seen from p)|) . (j27jl . that 

EZ n = (EZ n )(5) = - (\ogf n (e- 5 ))' , 5 > 0. (33) 

It is interesting to note that in the context of thermodynamics, the quan- 
tity \ogf n (e~ s ) has a meaning of the entropy of a system. This important 
fact that clarifies the choice of the free parameter was observed already by 
Khintchine ([14J, Chapter VI), in the course of his study of classic models of 
thermodynamics. 

From this point on, our study will be restricted to the three above mentioned 
classic combinatorial structures: multisets (weighted partitions), selections and 
assemblies. Recalling the forms of their generating functions i — 1,2,3 
(see [2]) and denoting ^(5) = /W (V* 5 ) , 5 > 0, i = 1, 2, 3, we obtain 



^)=n(i-e- 



k5\-b k 



k>l 



^ (2) w=n(i+e 



~kS\b k 



k>l 



^ 3 )(5)=exp($> e - fe5 ]. (34) 

yk>l 



Now it is easy to derive from (1281) and (|34|) that the following three types 
of distributions for the random variables ^Yj, in (|28j) : Negative Binomial 

(bk',e' 5k ^j , Binomial (bk, i+J-tk ) and Poisson (&/ c e- <5A: ), produce respectively 

cW, z = 1, 2, 3 in the representation f[3~Tj) . 

The representation (13T!) for assemblies was obtained in [9], while the one for 
general multisets and selections was obtained in [10] . The corresponding trun- 
cated generating functions ffl(z) are 
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/i i) w=n(i 



k=l 



z) b \ 



k=l 



=exp(£& fc s*J. (35) 



Consequently, in the three cases considered the equation (1321) takes the forms 
([IID-flSED derived from ([33]) 



» kb k e~ k ^ 

2^- ZIo> =n ' ( 36 ) 



k=l 



1 - e- k5: 



7(3)7 = ^ ( 3T ) 



fc=l 



1 + e" 5 " fc 



Y,kb k e- 5 " )k = n. (38) 
fe=i 



5 Preliminary asymptotic results 

In this section we find asymptotics for solutions to (|36l) - (!38l) . 

Lemma 2 Suppose that the sequence b k > 0, k > 1 is such that the asso- 
ciated Dirichlet generating function D satisfies the conditions (i) and (ii) of 
Theorem^ Then 

(i) As5-> 0+ 

^ (1) (5) =exp (AT(r)C(r + l)5~ r - D(0) log 5 + D'{0) + O (5 Co )), (39) 

f i2) {5) =exp (Ar(r)(l - 2~ r )C(r + l)cT r + D(0) log 2 + O (S Co )), (40) 

J^ (3) ((5) =exp (Ar(r)<T r + D(0) + O ($ C °)), (41) 

whereas asymptotic expressions for the derivatives 

^log^W(5)J , i = 1,2,3; k = 1,2,3 
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are given by the formal differentiation of the logarithms of tfJU\) -frfI\): 



(\ W 
log^ (1) (5)j = (-l) fc Ar(r + k)((r + l)5- r - k + (-!)*(* - l)D(0)r fc + O (8 c °- k ) , 



log^ (2) (<5) j = (-l) k AT(r + k)(l- 2- r )((r + l)<T r - fc + O (<5 C °" 
log^ (3) (5))) = (-l) k AT{r + k)S- r - k + O (8 c °~ k ) . (42) 



(ii) Each of the equations ([2f-(E2f has a unique solution <jW s^c/i £/ja£ 
<JW _>o, n^oo, i = 1,2,3. 

Moreover, 

(Hi) As n — > oo, 



<5« = (AT(r + l)C(r + l))^ 1 + ^U- 1 + OtrT 1 ^), 

r + 1 

u>/iere 



/?=<^' z/r " Co (43) 
^q-j-, otherwise] 



5i 2) = (Ar(r + l)(l-2- r )C(r + l)) r + 1 n-^T+0(n- 1 ~ /3 ) 



where /? = —5-; (44) 
r + 1 



<^ 3) = (AT(r + l))^n~^r +0 (V 1 ^) , 



where (3 = — . (45) 
r + 1 



fry) yls n — * oo, /^(e <5 ™ ) ), % = 1,2,3 have the asymptotic expansions of the 
right hand sides of / fffPj) - fiJUfy respectively, with 8 = 8$. 

Proof (i) First consider the case of weighted partitions. Following the Meinar- 
dus approach, we will use the fact that e~ u , u > 0, is the Mellin transform of 
the Gamma function: 



rv+ioo 

/ u~ s T(s)ds, u>0,v>0. (46) 

J v—ioo 



| pv+ioo 

2wi Jv 



14 



Expanding log in (154"|) as 



log rV>(S) = - £ h log(l - e~ 5k ) = E - E he- Skj 

k>l j>l J k>l 

and substituting f|46l) with t> = 1 + r gives 

log = — / S- s T(s)C(s + l)D(s)ds. (47) 

Z7TZ Jl+r-ioo 

By Meinardus' condition (i), the function D has a simple pole at r > with 
residue A, which says that the integrand in (14 7p has a simple pole at s = r 
with residue A5~ r T(r)((r + 1). Next, from the Laurent expansions at s — 
of the Riemann Zeta function £(s + 1) = - + 7 + . . . and the Gamma 
function T(s) = — 7+ . . . , where 7 is Euler's constant, and the Taylor series 
expansions at s = of the two remaining factors of the integrand in (|47[) . one 
concludes that the integrand has also a pole of a second order at s = with 
residue D'(0) — D(0)log5. We also recall that the only poles of T(s) are at 
s = —k, k = 0, 1, . . . Hence, in the complex domain —Cq < 3?(s) < 1 + r, 
with < Co < 1, the integrand has only two poles at and r with the above 
residuals. We now apply the residue theorem for the integrand in (147p . over 
the above domain. The assumption (jHJ) and the following two properties of 
Zeta and Gamma functions 

C(a + 1 + it) = O (\t\ C2 ) , t -> 00, C 2 >0, 



r(a +it) = O ^|t| Cs exp(-^|t|)^ , t -»• oo, C 3 > 0, 

uniformly in a, allow us to conclude that the integral of the integrand con- 
sidered, over the horizontal contour —Cq < 3?(s) < 1 + r, S(s) = t, tends to 
zero, as t — > 00, for any fixed 5. Thus, we are able to rewrite (T4Tj) as 



log f( 1 »((5) = AtT(r)((r+ 1) - D(0)\og5 + D'(0) 

1 r—Co+ioc 

+ — 6- s T(s)as + l)D(s)ds. (48) 

27T? J-Co-ioo 

Moreover, the previous two bounds and the bound (j3J) in Meinardus' condition 
(ii) imply that the integral in (HHl) is bounded by 



1 e— Co+ioo 

— / rT(a)C(a+l)I?(a)ds 



O (5 C ° l°° exp |t| Cl+C2+c ^^ 
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This proves (|39|). 



To prove the asymptotic formula for the first derivative I log T^ 1 ' (5) I , one 



has to differentiate (1481) with respect to 5 and then to estimate the resulting 
integral in the same way as above. Subsequent differentiations produce the 

asymptotic formulae for logjFW^) ) , k = 2,3. 



The proof of part (i) of the theorem for selections and assemblies is done in 
a similar way we now briefly describe. Following f[3~4"|) . the representation (|47l) 
conforms to 

1 fl+r+ioo 

log J™{S) = — / <TT( S )(1 - 2- s )C(s + l)D(s)ds (49) 

and 

log J* 3 \6) = — / S- s T(s)D(s)ds, (50) 

for all 5 > 0. Accordingly, the integrand in (1491) has a simple pole at s = r > 
with residue A5' r T(r)(l — 2~ r )((r + l), and a simple pole at s = with residue 
D(0) log 2, while the integrand in (150]) has two simple poles at s = r > and 
s = with residues v4<5~ r r(r) and -D(O) respectively. As a result, we obtain 
(J40D and dHJ. 

(ii) We see that the left hand sides of the equations (I36H38I) are decreasing as 
5 > in such a way that for a fixed n, in all the three cases the left hand 
sides tend to as 5 — > +oo, while as 5 — > the left hand sides tend to +oo, 
\ J2k=i and J2k=i respectively. We now make use of Theorem [2] to get 
a lower bound (1BT]) below on the sum J2k=i when the sequence obeys 
Meinardus' conditions (i) and (ii). We set a& := k~ r+1 hk, k > 1 and let -D(s) 
denote the Dirichlet series D(s) = J2k>i a kk~ s ■ Since C*o,r > 0, the function 
D satisfies the conditions of Wiener-Ikehara theorem, with the constant A as 
in Meinardus' condition (i). Consequently, 

n n kbk 

k=l k=l K 

from which it follows that for sufficiently large n, 

n 

E kb k > Bn, (51) 

k=l 

for some B > 1. This can be easily seen from the bound 



n 



A^i h.r — t-^i t.r Tr ^ K 

k=l K k=l K ^ k=L 
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with L < n such that U A > 1). Moreover, (fTTj) implies that the series 
linin^oo EZ^, % = 1, 2, 3 converge for any positive 5. 

Combining the above facts, we conclude that each of the equations (I36H38I) has 
a unique solution for sufficiently large n and that the solutions 6$ — > 0, n — > 
oo, % = 1, 2, 3. 



(iii) We firstly show that in all three cases, 



EZ« = (-log^«(e" 



5=SS ) 



+ e ( %), e<*>(n) 



n^oo, i = 1, 2, 3. (52) 



r + 2 



In the case of weighted partitions, setting 5 n = n 2 ( r + 1 ) gives for sufficiently 
large n 



°° hhi p~ k5 n ( CO \ I CO 

E = °[ E =0 E ^ r+1 e- M " 

O ( J°°^ x r+l e- xK dx^j -»• 0, n -> oo, (53) 

where we have employed fill I) and the fact that n<5 n — > oo, n — > oo. From 
fj42l with k — 1 we deduce that that 5^ > <5 n for large enough n, which 
implies that (1531) is valid with <5 n replaced by 5^>. This proves fl52|) for the 
case considered. Consequently, the equation ( l36l) can be rewritten as 



AT(r + l)C(r + l)^)"^ 1 + IW?) -1 + O ((^ 1) ) c<) - 1 ) + e«(n) = n, 

e (1) (^)^0, n^oo. (54) 

We outline here the method of solution for asymptotic equations of the form 
(J|)J) common in applications of Khintchine's method. Denoting the constant 
coefficient h := AY{r + l)£(r + 1), (jSJ) implies 

* + ^X^f + O + o ((6^) = n (^)) r+1 . (55) 

Since 8$ — > 0, n — > oo, we obtain from fl55l) that #W ~ h~n~~ , n — > oo. 
Based on this fact and the fact that < C < 1, we get 



h r + 1 n r + 1 + 



+ 



r + 1 
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where (3 is as in (H3|) . For selections and assemblies the analogs of (I55I) will be 
respectively 



^o((^) C ° +r ) + o((^) r+1 )=.(^)) 
h = AT(r + 1)(1 - 2" r )C(r + 1) 



r+l 



and 



ft + o((^) G " + ') + o((^)' +I )= n (^)' +1 , 
h = AT(r + 1). 

Now the same reasoning as for weighted partitions leads to the solutions f|44|) . 

<m>- 

(iv) In the case of weighted partitions, we have 

log/«(e- 5 " ) ) = logJ- 1 (5i 1 ))+ £ b k log(l-e- k5 ^), 

k>n+l 

where 

O ( E &*e~ fe *"' ] = o (n r+1 exp (-n^)) 



E Wog(i-^"') 

fc>n+l 



v fc>n+l 

— > OO, 



by the argument giving ( l53l) . The proof of the remaining parts of the assertion 
(iv) is similar. ■ 

Remark As we mentioned before, Meinardus' proof (see pQ) of Theorem [1] 
relies on application of the saddle point method. In accordance with the prin- 
ciple of the method, the value in question is expressed as 

c « = _L f 1/2 ^( 1 )(r)e n5+2 ™ a cia, r = 6 + 2ma, (56) 
2m J -i/2 

by virtue of the Cauchy integral theorem. Here the free parameter 6 is chosen 
as the minimal value of the function exp (AT(r)((r + l)5~ r + nS) viewed as 
an approximation of the absolute value of the integrand in ( J56l) at a = 0. This 
gives 5 = h~n~~~ , h = AT(r + 1)C(?" + 1) which is the principal term of 
the solution S^' of fl36l) . It can be seen that, stemming from this choice of 
the free parameter, the subsequent steps of Meinardus' proof are considerably 
more complicated compared with ours. Also note that our choice of the free 
parameter is in the core of our ability to weaken Meinardus' condition (Hi). 



Our next assertion reveals that the function 9?(g(r)) — g(5) in the left hand 
side of the Meinardus' condition (Hi) is inherent in the employed probabilistic 
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method: the function provides an upper bound for the rate of exponential 
decay of the absolute value of the characteristic function n in (I29I) . as n — > oo, 
for all three types of random structures considered. The bounds obtained in 
the forthcoming lemma are used in the the proof of our local limit theorem, 
Theorem HI 

Recall that the function g{r) is defined by (j3J) and (jSJ). 
Lemma 3 Denote 



V(a) = V(a; 5) = ^(g(r)) - g(5), r = 5 + 2nia 

oo 

= -2 he~ kS sin 2 (neck), 5 > 0, aGl 
fc=i 

and let 5 = 8%\ i = 1,2,3 be the unique solutions of the equations (f375j ) -(fff5 P 
respectively. Then, for all aGl, 

I0i i) («)l<(l + ^ ) )exp(^^), eS? = c 0( a )_> O , n - oo, » = 1, 2, 3, 

where V^\a) = V{pt]8$) and the constants M^, i = 1,2,3 are as in condi- 
tion (iii f ) of Theorem^ 

Proof From (I29I) we have for n > 1 and 5 > fixed, 

|0«(a)| = exp ^ (log/«(e-)) - log^e" 5 )) := e v » ^ 

r = 5 + 2mia. (57) 

Using fl35|) we now find bounds for V^^(a;5) expressed via V^ l '(a;8) in the 
three cases considered. By the definition of V^(a; 5) as given in fl57|) . we have 
for weighted partitions, 



V«(a;*) = R(-2>k>g( T 



1 - e 



-rfc ' 



fe=l 
1 ™ 

log 

z fc=i 



e 

,-fe<5 



1 -2e- kS cos(2vraA;) + e- 



(1 - e- Sk ) 2 

1 » / 4e- 5fc sin 2 (7raA;) 

> Ofc lOg 1 H ; 

2 ^ 6 V 1 - 



k=l 
1 n 



-Sk\2 



< — & fc log ( 1 + 4e _<5fc sin 2 (vraA;) 
^ fe=i ^ 

lnp- 5 n 

<--|-^6 fc e- 5fc sin 2 (7raA;), 5 > 0, a£ 
2 k=i 
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where the last inequality is due to the fact that log(l + x) > (^f^)a;, < 
x < 4. 

For selections, in a similar manner, 

^)(a;tf) = »(g6 Jk tog(i±^)) 

1 » /l + 2e-^cos(2vraA;) + e~ 2fe5 \ 

-£ 6tlog ( aTF^p j 

1 » / 4e- fc5 sin 2 (7raA;)\ 

= 2g 1 6 * log l 1 " (1 + e-^ J 

_1 " 4e- fc<5 sin 2 (7ra£;) 
"~2^i * (1 + e^) 2 

1 n 

< --Y.he^ 5 sin 2 (nak), 5 > 0, «6l 

2 fe=i 

Here the first inequality is due to the fact that — log(l — x) > x, < x < 1. 
For assemblies, we get straightforwardly 

n 

V^ 3 \a;S) = -2]T& fc e~ fe<5 sin 2 (7ra£;), 5 > 0, a £ R. 
k=i 

Finally, setting 5 = 8$ in the above three expressions, the argument resulting 
in ( 153]) implies that in all three cases, 

K W («;^)<^^ + ^, 6« = e«(«)^0, n^oc, (58) 
uniformly for all aeK. This completes the proof. ■ 



6 The local limit theorem and completion of the proof 

Local limit theorems are viewed as the main ingredient of the Khintchine 
method. Theorem [4] below says that a local limit theorem holds for all three 
types of structures obeying the conditions of our Theorem [31 

Theorem 4 (Local limit theorem) . Let 8$, i = 1,2,3 denote the solu- 
tions to the equations fifty) -( fifty) respectively and let the random variables 
n > 1 be defined as in fifty) , where the random variables have distri- 
butions given in the paragraph following fifffy . Assume that condition (Hi') of 
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Theorem^ holds for i = 1, 2, 3. Then, 

nZ$=n)~ , 1 ^^_(^)) 1+r/2 ; * = 1,2,3, 

27rVar(ZL l) ) V 27ri V 



2 

itrai/i constants defined by 

^ 1) = Ar(r + 2)C(r + l), 
f^ 2) = A(l - 2- r )T(r + 2)C(r + 1) 

and 

Kf ] = AY{r + 2). 

Proof We will find asymptotics for F(Z n = n) as n — > oo for the three types 
of random structures. Following the pattern of the Khintchine method (see 
e.g. [9lfT0] ). we set 5 = 5 n in (l28l) and (l29~I) and define «o = «o(^) to be 



r+2 



a = ( 5F TTT log 2 n. (59) 
Then we have 



F(Z n = n)= f 1/2 <f) n (®)e- 2nma da = h + J 2 , (60) 

J -1/2 

where I\ = I\{n) and J 2 = I<i{n) are defined to be 

I x = r° <f) n (a)e- 2mna da (61) 

J—ao 

and 

n (a)e" 27rma da + / n (a) e - 2 ™ a cia. (62) 

-1/2 ./ceo 

Defining 5„ and T n by 

^ = ^ (log/n(e- 5 ")) , T n = (log /n (e- 5 ")) , (63) 

for n fixed we have the expansion 



n (a)e" 2 ™ a = exp [2mot{EZ n - n) - 2ix 2 u 2 B 2 n + 0(a 3 )T n 

= exp (-27r 2 a 2 J B 2 + 0(a 3 )T n ) , a -> 0. (64) 

It can be checked that B 2 n = VarZ n and T n = £? =1 E(^ - EY}) 3 , by the 
argument leading to ( 1331) . 
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Now (1421) in Lemma [2] and (163]) tell us that for all structures considered 

(Blf ] ~ K%\5^yr-\ (65) 

and 



where > 0, i = 1,2,3, are constants depending on the type of the 

structure and while i — 1, 2, 3 are as in the statement of the theorem. 

Therefore, considering (1591) . we find that in all three cases, 

lim Bic? n = oo and lim T n af, = 0. (66) 

n— »oo u n^oo u 

Combining (jM|) with (1661) . we deduce that 

0„(a)e" 2 ™ a ~ exp (-2tt VB^) , n -> oo, \a\ < a . (67) 
Finally, using (IBTI) . (|66|) and (1671) gives us 



exp {-2n 2 a 2 B 2 n )da ~ (27TEJ" 1 / ( 

-cm J —oo 



'da = , n — » oo. 

(68) 

The next step of the proof is to show that J 2 = o(Ji), n — > 00. At this step 
condition (m') of Theorem [3] plays a key role. Because of the asymptotic 

\l 8$ = o(a^), n — > 00, i = 1, 2, 3, we can use condition (m') to bound the 
quantity V^'(a) defined in Lemma [3] by 

V {€} {a) < (\+ r - + ^jM®\og8%\ a < \a\ < 1/2, 1 = 1,2,3. 

Hence, Lemma [3] and the fact that in condition (iii'), e > give 



^(a)| = o((5«) 1+i )(l + 6«), a < \a\ < 1/2, n - 00, * = 1,2,3. 



From the definition (1621) and the asymptotic (1651) we have 



h = o(m 1+r/2 )=o(I 1 ). (69) 



Lastly from (160|) . (1651) . (1691) and (1681) . we derive the following asymptotic ex- 
pression for F(Zj?' = n), i = 1, 2, 3: 



l+r/2 



, n — > 00. 
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To complete the proof of Theorem [3] it is left to substitute the asymptotic 
expressions implied by our results for the three factors in the representation 
(1311) when 5 = <jW. 



7 Concluding remarks 

(i) Under the stated conditions on parameters bf., the asymptotic formulae 
©, (l20|) - fl2~T|) . for c^, % — 1, 2, 3 have a striking similarity, all of them being 
of the form: 



where we have denoted by xi, X2, X3 the constants that depend on the type of 
a structure and it parameters. A simple analysis that takes into account that 
C(r + 1) > 1 reveals that, asymptotically in n, c£> > and > where 
the first fact follows obviously from the definition of selections. 

(ii) The following observation is also in order. It turns out that each one of 
the three combinatorial structures obeying the above conditions behave very 
much alike the one with parameters b^ = /c r_1 , r > 0. According to the 
classification suggested in [I] the latter structures are called expansive. In 
this respect, combinatorial structures obeying Meinardus' conditions (as well 
as their extensions as defined in the present paper) can be viewed as quasi- 
expansive. 

(iii) We hope that the approach of this paper can be applied as well to other 
enumeration problems, in particular to enumeration of structures with con- 
straints on the number of summands (components) (see e.g. [19J). 
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